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Abstract 



The production of forward jets of transverse energy Et — Q and large momentum 
fraction Xj e t 3> x is calculated in next-to-leading order including consistently direct 
and resolved virtual photon contributions. The predictions are compared to recent 
ZEUS and HI data. Good agreement with the data is found. 



1 Introduction 



The cross section for forward jet production in deep inelastic scattering (DIS) has been 
proposed as a particularly sensitive means to investigate the parton dynamics at small 
x [|]]. Analytic calculations based on the BFKL equation in the leading-logarithmic 
approximation show a strong rise of this cross section with decreasing x [|[ and were found 
in reasonable agreement Q with the first data from the HI collaboration at HERA ||. 
More recent measurements of the forward cross section, based on an order of magnitude 
increased statistics compared to [f| , have been presented recently by the ZEUS @ and the 
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HI [0 collaborations confirming the earlier findings ||. Monte Carlo generators based 
on direct photon interactions (DIR) calculated from leading order (LO) 0(a s ) matrix 
elements together with leading-logarithm 

parton showers disagree with the measured jet cross section || |7j by an appreciable factor. 
Also next-to- leading order (NLO), i.e. 0(a 2 ), calculations predict too small forward cross 
sections at small x as already shown by Mirkes and Zeppenfeld || using their MEPJET 
program || when comparing to the data in ||. This has been confirmed also with the 
new data in P, H . 

A similar deficiency between NLO calculations and measured data occurs for the dijet 
rate in the region E\ > Q 2 [|10J, where Et is the transverse energy of the produced jets 
and Q 2 is the usual squared lepton momentum transfer. This kinematic range is also 
relevant for the forward jet production, as will be seen later. The region of small enough 
Q 2 is the photoproduction regime where the virtual photon resolves into partons. Indeed, 
introducing a resolved photon contribution, the measured dijet rate and the forward jet 
cross section can be described satisfactorily |]TTJ, [L0|, |7j concerning the shape of the cross 



section as a function of x as well as the absolute normalization. This description is based 



on the Monte Carlo program RAPGAP [12| which includes a resolved photon contribution 
in addition to the direct process, which both are evaluated with LO matrix elements with 
additional emissions in the initial and final state generated by parton showers together 
with subsequent hadronization. 

The dijet rate has been calculated also in NLO including direct and resolved photon 
contributions In order to avoid double counting in the full NLO calculation, the 

contribution from the virtual photon splitting into qq pairs, where either the quark or the 
antiquark subsequently interacts with a parton originating from the proton, had to be 



subtracted JR], similar as is done in the NLO theory for the photoproduction of jets |T5 



The subtracted terms in the NLO direct contribution are part of the parton distribution 
functions (PDF's) of the virtual photon and appear in the resolved contribution in an 
evolved form. With this procedure the whole cross section for two-jet production, which 
is a superposition of the direct and resolved contributions minus the photon splitting 
piece, becomes to a large extent independent of the factorization scale at the photon 
vertex. This full NLO calculation of the dijet rate agreed well with the HI data over the 
full Q 2 domain, 5 < Q 2 < 100 GeV 2 , and the x domain, 10~ 4 < x < 10~ 2 , and for jet 
transverse momenta E\ > Q 2 |ff| |13|. 

In this work we want to present the results of a calculation of the forward jet cross section 
on the basis of the NLO theory used for the dijet rate. Although the kinematic constraints, 
very low x and E\jQ 2 of order one, are rather similar, it is not obvious that the calculated 
cross sections will agree with the recent ZEUS || and HI J7| experimental results. 

After some comparisons with the MEPJET results [0] to make sure that our DIS jet program, 



called JetViP ||13|| , gives the same results under identical kinematical conditions we shall 
give our results with the experimental cuts of the ZEUS M and HI |7J] analysis. We close 
with a short summary and an outlook to future studies. 
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2 Comparisons and Results 



2.1 Comparison with MEPJET 

Before we present our results with the ZEUS and HI kinematical constraints for select- 
ing the forward jets we performed a check of our program JetViP with the forward jet 
kinematics by comparing with the NLO results of Mirkes and Zeppenfeld ||, who have 
produced their results with the fixed order program MEPJET ||, which only includes direct 
photon contributions. We have chosen the same kinematical cuts, which differ somewhat 
from the cuts used in the ZEUS [[| and HI Q analyses. 

The 0(a s ) results are obtained taking the Gliick, Reya and Vogt (GRV) LO proton PDF's 
T6fl together with the one-loop formula for a s . For the 0(a 2 ) results we employ the GRV 



higher order PDF's together with the two-loop running a s formula. We take Nf = 5 and 
match the strong coupling at the charm and bottom thresholds fiR = m c , mt,, respectively. 

Jets are defined in the laboratory frame using the cone algorithm with the opening angle 



(A77) 2 + (A0) 2 = AR < 1 in the so-called E-scheme. In this scheme the four-vector 
of the combined jet is given as the sum of the four-vectors of the two partons. The 
differences of pseudorapidities and azimuthal angles with respect to the jet direction are 
Arj and A(j). All jets have to fulfill \r]\ < 3.5 and Et, E^ > 4 GeV, where the index B refers 
to quantities in the Breit frame, rj and Et are measured in the HERA laboratory frame. 
Additional cuts are made for events which contain a forward jet. These requirements are 
1.735 < r] < 2.90 and E T > 5 GeV with 

p z /E P > 0.05 , 0.5 < E T /Q 2 < 4 . (1) 

The x variable is restricted to the small-re region of x < 0.004. The cuts on the electron 
variables are Q 2 > 8 GeV 2 , y > 0.1, E' > 11 GeV and 0' e G [160°, 173.5°]. The electron 
and proton energies are E e = 27.5 GeV and E P = 820 GeV, respectively. The positive 
z-direction is the direction of the incoming proton momentum. 

The renormalization (/xr) and factorization scales (hf) are taken equal and are identi- 
fied with the sum = /ip = where fc|?(z) and p?, the parton's transverse 
momentum, are related in the Breit frame by 

[k^)f = 2^ 2 (1 - cos**) = - 2 [^(*)] 2 , (2) 

1 + COS (7 ip 

where 9 ip is the angle between the parton and the proton direction in the Breit system. 
In LO one-jet production, i.e., in the naive parton model limit, k^(i) = Q. With these 
constraints we obtain the cross sections in Tab. 1, where our JetViP results are compared 
to results from ||, referenced as MEPJET in the table. In addition to the forward jet 
cross sections we also list the full 2-jet, inclusive 2-jet and exclusive 3-jet cross sections. 
Compared to these cross sections, the kinematic constraints defining the forward cross 
section lead to considerable reductions. Furthermore we notice that the NLO corrections 
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Table 1: Jet cross sections in the forward region compared to MEPJET results. 



Contribution 


MEPJET 


JetViP 


relat. difference 


0(a.) 2 jet 


2120 pb 


2203 pb 


+4% 


same + forward jet 


18.9 pb 


20.0 pb 


+6 % 


0(a 2 s ) 2 jet inclusive 


2400 pb 


2371 pb 


-1% 


same + forward jet 


83.8 pb 


89.0 pb 


+6% 


0(0*1) 3 jet exclusive 


210 pb 


207 pb 


-1% 


same + forward jet 


14.8 pb 


14.5 pb 


-2% 



to the forward jet cross sections are large in agreement with ||. The 0(oP s ) single-jet 
cross section is not considered since it vanishes if a forward jet is required, due to the 
kinematical restrictions of the phase space. 

The numbers for MEPJET are taken from ref. ||. They differ by a few percent from our 
JetViP results. This is due to a different implementation. In JetViP the azimuthal (0) 
dependence of the jet with respect to the electron plane is integrated out in the hadronic 
center-of-mass system, whereas in MEPJET this <fi dependence is included in this frame and 
then integrated in the HERA laboratory system. These terms which originate from the 
interference of the longitudinal and transverse virtual photon polarization (~ cos </>) and 
from the transverse linear photon polarization (~ cos 20) vanish for Q 2 — > ||17|| . Since in 



our case the virtuality Q 2 is not very large, the contribution of the azimuthal dependent 
terms is small, which leads to small deviations between MEPJET and JetViP results. In 
addition, when integrating over the full phase space, it does not matter in which system 
the 4> integration is performed, so that the observed difference is essentially due to the 
phase space restrictions in the forward jet selection. 

Next we have evaluated the forward jet cross section including a resolved virtual photon 
contribution. Since the Q 2 is fairly large one might think that introducing a resolved 
virtual photon component is not necessary, because it is equivalent to a contribution 
generated in the NLO correction to the direct cross section. However, the NLO resolved 
cross section introduces additional higher order terms that are not contained in the NLO 
direct cross section, as we will see below. The kinematical cuts for selecting the forward 
jet are the same as for the comparison with MEPJET. The resolved cross section in LO and 
up to NLO is calculated as described in our earlier work jy| and which is incorporated in 
the JetViP program [13]. The PDF's of the virtual photon are taken from ||18|| , specifically 



we took the version SaSID, which was transformed to the MS scheme (see |13J). For the 



LO resolved cross section it would be more appropriate to choose SaSID without the 
transformation to the MS scheme. This would increase the total LO cross section in Tab. 



2 by 10%. As in |13j we subtracted a term originating from the 7* — > qq splitting in 



4 



Table 2: Resolved component in the forward region. 



Contribution 


0(a s ) 2 jet 


0(a 2 s ) 2 jet incl. 


Resolved 


61.3 ±0.5 pb 


109.6 ±0.7 pb 


7* — > qq splitting 




-52.9 ±0.3 pb 


DIS Direct 


18.7 ±0.2 pb 


89.6 ±0.3 pb 


Sum 


80.0 ±0.5 pb 


146.3 ± 0.8 pb 



the NLO direct photon matrix elements in order to avoid double counting at the NLO 
level. The SaS parametrizations of the virtual photon PDF vanish if the virtuality Q 2 is 
larger than the factorization scale pip. Therefore we have chosen fx 2 R = fx 2 F = Q 2 + (E F ) 2 . 
This enforces the virtual photon to be present since now always Q 2 < fj, 2 F . The results 
for the various components of the forward cross section are summarized in Tab. 2. From 
these results we observe the following. First, the sum of the LO direct and resolved 
contributions coincide within 10% with the NLO direct cross section. Second, adding the 
subtracted NLO direct contribution, which is the NLO direct minus the contribution from 
the photon splitting term (given with the minus sign in Tab. 2) to the NLO resolved cross 
section, leads to a large correction of about 60% if compared to the full NLO direct cross 
section. This increase has two sources. First, the LO resolved cross section is 15% larger 
than the subtracted photon splitting term. Part of this increase is due to the evolution 
of the PDF's of the photon. The other part comes from the gluon component of the 
photon PDF, which amounts to 7.3 pb. Second, the NLO corrections to the resolved 
cross section give a further increase as compared to the LO result of approximately 80%, 
which originates from the NLO corrections to the resolved matrix elements. 

We conclude, that the NLO resolved contribution supplies higher order terms in two ways, 
first through the NLO corrections in the hard scattering cross section and second in the 
leading logarithmic approximation by evolving the PDF's of the virtual photon to the 
chosen factorization scale. This way we sum the logarithms in E T /Q 2 , which, however, 
in the considered kinematical region is not an important effect numerically, as we have 
seen in Tab. 2. Therefore, the enhancement of the NLO direct cross section through 
inclusion of resolved processes in NLO is mainly due to the convolution of the point- 
like term in the photon PDF with the NLO resolved matrix elements, which gives an 
approximation to the NNLO direct cross section without resolved contributions. One of 
the dominant contributions to the forward cross section is shown in Fig. 1 (left part), where 
the photon splitting term is convoluted with a matrix element that provides two gluons 
in the final state. This way one gluon rung is added to the gluon ladder as compared to 
the corresponding NLO direct cross section, shown on the right of Fig. 1. The additional 
gluon in the NLO resolved term is in our approach calculated from perturbative QCD, 
producing an additional term which makes a contribution to the forward jet. In the 
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Figure 1: Diagrams contributing to resolved (left) and direct (right) processes in the 
forward region. 

BFKL approach for the forward jet cross section || || this extra gluon is part of the 
BFKL evolution. In contrast, the NLO direct term in Fig. 1 (right part) contains also 
additional gluons in the DGLAP evolution of the proton PDF, which, however, are not 
resolved, i.e., go to the proton remnant. Another way to generate a larger forward jet 
cross section would be to go to the NNLO corrections of the direct cross section, which 
has not been done yet. We think that including the NLO resolved component produces a 
reasonable approximation to this NNLO cross section. Such a correspondence is present 
at one order lower. As already remarked above, the superposition of the LO direct and 
resolved cross section is almost equal to the NLO direct cross section. 



2.2 Comparison with ZEUS and HI Data 

For the comparison with the 1995 ZEUS || and the 1994 HI |7J forward jet cross section 
data we calculated the NLO cross section with slightly different input and in particular 
with the exact kinematical constraints for the forward jet selection as used in the two 
experiments. 



As proton PDF's we apply now the CTEQ4M parametrization |19| with the two-loop a s . 
We take Nf = 5 as before and match the value of a s at the thresholds /ir = m c , with 
a Ajjg as used in CTEQ4M. Jets are defined with the cone algorithm in the HERA frame 
as described above, except that the axis of the jet is calculated now as the transverse 
energy weighted mean of rj and <fi of the two partons or jets belonging to the combined 
jet. This kind of jet definition was also applied in the experimental jet analysis. As scales 
we choose /i 2 = [i\ = \i\ = M 2 + Q 2 with a fixed M 2 = 50 GeV 2 related to the mean 
E T of the forward jet. We take this fixed value of M instead of Et for technical reasons, 
since the calculations in JetViP start from the hadronic c.m.s.. The choice n 2 F > Q 2 is 
mandatory if we want to include a resolved contribution. Another choice of scale would 
be [If = Et or /ip = M. In the case > 1 only for E T /Q 2 > 1, which covers only 

part of the ZEUS kinematical range. To have a resolved cross section in all E\jQ 2 bins 
we consider the choice n 2 F = M 2 + Q 2 more appropriate. 

In the two experiments the forward jet selection criteria are different. In the ZEUS 
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Figure 2: Dijet cross section in the forward region compared to ZEUS data, (a) NLO 
DIS, E T > 5 GeV (b) NLO DIR 5 +RES, E T > 5 GeV. 



experiment the kinematical constraints are: E' e > 10 GeV, y > 0.1, r\ > 2.6 and Et > 
5 GeV with the forward jet constraints Xj et = Ej et /Ep > 0.036, 0.5 < E\jQ 2 < 2, 
Pzjet > an d 4.5 x 10~ 4 < x < 4.5 x 10 -2 . Our results for the forward jet cross section 
under these ZEUS kinematical conditions are shown in Fig. 2 a,b. In Fig. 2 a we plotted 
the full 0(a 2 ) inclusive two-jet cross section (DIR) as a function of x for three different 
scales fi 2 = 3M 2 + Q 2 , M 2 + Q 2 and M 2 /3 + Q 2 and compared them with the measured 
points from ZEUS 0. As to be expected the calculated NLO direct cross section is by a 
factor 2 to 4 too small compared to the data. The variation inside the assumed range of 
scales is small, so that also with a reasonable change of scales we can not get agreement 
with the data. In Fig. 2 b we show the corresponding forward jet cross sections with 
the NLO resolved contribution included, as described in the previous subsection, again 
for the three different scales \x as in Fig. 2 a. Now we find good agreement with the 
ZEUS data. The scale variation of the calculated cross section is larger than in Fig. 2 
a. In particular, the largest scale gives now the largest cross section opposite to what 
is observed in Fig. 2 a. This different scale variation comes primarily from the scale 
dependence of the virtual photon PDF. This factorization scale variation is supposed to 
be compensated between the LO resolved and the NLO direct contribution but not for 
the NLO resolved contribution. This could only occur if we could include the NNLO 
direct contribution which, however, is not available. Since the NLO resolved contribution 
for the forward jet is rather large, as discussed above, this scale dependence can not be 
avoided. On the other hand, the scale dependence is not so large that we must fear our 
results not to be trustworthy. In Fig. 2 b the cross section is labeled DIR5+RES, where 
DIR5 stands for NLO direct minus the photon-quark-antiquark splitting term. We have 
calculated the forward jet cross section also with the scale [i 2 F = M 2 . For this choice we 
obtain a somewhat smaller cross section which is approximately equal to the cross section 
with the scale fi% = M 2 /3 + Q 2 . 

In the HI experiment the forward jets are selected with the kinematical cuts E' e > 
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Figure 3: Dijet cross section in the forward region compared to HI data, (a) NLO DIS, 
E T > 3.5 GeV (b) NLO DLR S +RES, E T > 3.5 GeV; (c) NLO DIS, E T > 5.0 GeV (d) 
NLO DIR5+RES, E T > 5.0 GeV. 



11 GeV, y > 0.1, 160° < 9' e < 173°, 1.735 < rj < 2.794 (this corresponds to 7° < 
9 jet < 20°), E T > 3.5 (5.0) GeV, x jet > 0.035 and 0.5 < E T /Q 2 < 2. This corresponds 
approximately to the Q 2 range 5 < Q 2 < 100 GeV 2 . The forward jet cross section 
is measured for various x bins ranging from 1.0 x 10~ 4 to 4.0 x 10~ 3 . Otherwise the 
calculated forward cross section is obtained under the same assumptions as for the ZEUS 
selection cuts. In Fig. 3 a,b,c,d we show the results compared to the HI data obtained 
with two E T cuts in the HERA system, E T > 3.5 GeV (Fig. 3 a,b) and E T > 5.0 GeV 
(Fig. 3 c,d). In the plots on the left (Fig. 3 a,c) the data are compared with the pure NLO 
direct prediction, which turns out to be too small by a similar factor as observed in the 
comparison with the ZEUS data. In Fig. 3 b,c the forward jet cross section is plotted with 
the NLO resolved contribution included in the way described above. For both Et cuts, 
Et > 3.5 GeV (in Fig. 3b) and Et > 5.0 GeV (Fig. 3d), we find good agreement with the 
1994 HI data inside the scale variation window M 2 /3 + Q 2 < /i 2 < 3M 2 + Q 2 . Compared 
to the ZEUS data the HI measurements extend down to smaller x. In the lowest x bin, 
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1.0 x < x < 5.0 x 10~ 4 the forward jet cross section has a dip, which is also reproduced 
in the theoretical calculation and which is due to the kinematical constraints for selecting 
forward jets. 

We conclude that the NLO theory with a resolved virtual photon contribution gives a good 
description of both, the ZEUS and the HI forward jet data. It is important, that both 
components, the direct and the resolved one, are calculated up to NLO. A LO calculation 
of both components would fall short of the experimental data, as it is clear from the 
results presented in Tab. 2, where we compared LO and NLO results. 

We remark that the forward jet cross sections as measured by ZEUS and HI are obtained 
at the hadron level, i.e., the jets are constructed from measured hadron momenta using 
the same cone algorithm. In our NLO calculation the jets are combined from partons with 
the same jet algorithm. The size of the corrections from hadron to parton level has been 
studied by the ZEUS collaboration using several Monte Carlo simulation programs 
with the result that for the models which account very well for the ZEUS forward jet 
cross sections the correction factors are close to unity for all x values considered in the 
analysis. In the HI work similar results are reported 0. 

In the HI forward jet data are also successfully described with the RAPGAP model [12 



which includes a direct and a resolved component, both in LO, and with a similar scale 
fi as we have used. In addition this model has leading logarithm parton showers in the 
initial and final state built in. We think that these parton shower contributions produce 
the higher order effects which we found necessary to account for the correct normalization 
and the x dependence of the forward jet cross section data. 



3 Concluding Remarks 

We conclude that the measurements of the forward jet cross section presented recently 
by the ZEUS and HI collaborations can be described very well by the NLO theory with 
direct and resolved virtual photon contributions added in a consistent way. The theory 
shows good agreement with the data with respect to the normalization and also the 
functional dependence with decreasing x. Whereas this variation of the cross section with 
x is also compatible with the NLO predictions based on direct photons, for the correct 
normalization the resolved component up to NLO is needed. To avoid double counting 
the 7* — > qq splitting term is removed from the NLO direct contribution. 

In contrast, LO BFKL predictions || yield much larger forward cross sections than the 
data m 0. These calculations suffer, however, from several deficiencies. They are asymp- 
totic and do not contain the correct kinematic constraints of the produced jets ^ . 



Furthermore they do not allow the implementation of a jet algorithm as used in the exper- 
imental analysis. Also NLO ln(l/x) terms in the BFKL kernel []ST] predict large negative 



corrections which are expected to reduce the forward cross section as well. When all these 
points are taken into account the BFKL approach may give an equally well description 
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of the forward jet data. Even if this is the case, it is clear from this work that the BFKL 
theory is not the only theoretical approach that describes the forward jet cross sections. 
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